Simple description of metastable phase decay 
on wide spectrum of activities 

V.Kurasov 
Victor.Kurasov@pobox.spbu.ru 



One of the typical examples of the system evolution in the first order 
phase transition is the metastable phase decay in the closed system. After 
the instantaneous creation of initial supersaturation $ the external influence 
doesn't take place and the evolution of the supersaturation £ in the system 
occurs only due to the nucleation kinetics. The supercritical embryos of a 
new liquid phase, i.e. the droplets, appear mainly on heterogeneous centers 
which inevitably exist in any real system. The number of heterogeneous cen- 
ters slightly changes in time, but their activation barriers can attain rather 
arbitrary values. Then one can speak about the distribution r]t t of the het- 
erogeneous center numbers over the heights AF of activation barriers (i.e 
over activities w of heterogeneous centers). 

It is natural to suppose that the mentioned distribution is rather smooth 
as a function of w and covers some region of the activation barriers heights. 
The subject of current publication is to give the most simple variant for 
the nucleation description (the full theory is given in [|1]). The property of 
activity spectrum to be the wide one will be seen a few moments later. 

The nucleation description has to take into account the processes of the 
mother phase consumption and the processes of exhaustion of heterogeneous 
centers. Then one comes to the non-linear system of condensation equations 
where the behavior of ( depended on / in all previous moments of time. The 
approximate method to solve this non-linear problem is given below. 

We consider the homogeneous system of a unit volume. The regime of 
the droplets growth is supposed to be the free molecular one, the process will 
be an isothermal one. 
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One can prove analytically the validity of quasistationary approximation 
for the nucleation rate for all situations when the critical embryo has at least 
few dozens of molecules. 

Under these conditions one can also prove that the main role in vapor con- 
sumption belongs to the supercritical embryos (which are called as droplets). 

The rate of nucleation / is the intensity of the droplets formation. To 
characterize the droplet with v molecules it is convenient to introduce the 
linear size p = u 1 ^ 3 . The rate of p growth is (/t with some characteristic 
constant r for all sizes of droplets. 

One can easily extract the main dependencies in the nucleation rate be- 
havior over the supersaturation ( and over the free (unoccupied by droplets) 
number 7] of the heterogeneous centers: / ~ exp(— AF)i], where AF is the 
height of the activation barrier taken in thermal units. The end of nucleation 
corresponds to the relatively small fall of supersaturation from initial value 
$. Then one can linearize AF as a function of ( at the nucleation period 
and come to the following approximation 

I(C,V) = I($,Vtot))(v/Vtot) ex P (r(c - $)/<£>) 

where T = —^dAF/d(\(; = q>. One can also see that during the nucleation the 
rate of droplets growth (/t is approximately constant in time $/r. 

To get the behavior of the supersaturation one has to know C(£). It can 
be gotten by the balance equation 

C = $ - G/ noo 

where G is the number of the molecules in the liquid phase and is the 
number of the molecules in the saturated vapor. For G one can get G = 
f dwg(w), where the integral is taken over the whole spectrum of activities 
and g(w) is the number of molecules in the droplets on the centers with 
activity w. The last value is equal to 

g(w)= [\t-t') 3 (-) 3 I(t',w)dt' 

JO T 

Having mentioned that the main dependence of I over the activity w is 
~ exp(— w) one can come to 

G= fdw [\t-tfe W (-w+w*)^^(-f^ 

J JO Vtot{W*) T <P 
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with some parameter w*. In the last equation there exists i](w,t f ). It can be 
given by the following expression 

7](w,t) = r] tot {w) exp(- / exp(-w+w*) — r exp(r )dt) 

The last two equations form the closed system of condensation equations. 

One can see that these equations are essentially non-linear ones. 

The choice of w* is rather arbitrary. ^From the last equation one can see 
that if w* is chosen to satisfy condition r](w*,t = oo) = r] tot {w lt )/2 then all 
centers with w > w* + 1 remain practically unexhausted and all centers with 
w < w* — 1 are practically exhausted. So, the interesting region [w* — 1, 
is rather narrow. One can put in [w* — 1, w* + 1] the value i] tot (w) to ri tot (w*). 

Distribution r/ tot (w) has to be limited from below to ensure the finite value 
of integral nucleation rate. Denote by Wq the minimum of w. If w* — Wq > 
(3-^-4) then one can say that the relatively "wide spectrum" of activities takes 
place in nucleation. In the opposite situation one can consider nucleation as 
taking place on the similar heterogeneous centers. 

Consider the situation of wide spectrum. We need to know the behavior 
of supersaturation at the end of nucleation when the fall of supersaturation 
stops the nucleation. All centers with w > w + = + 2 will be exhausted 
until the first quarter of the nucleation period duration. The droplets formed 
on these centers can be regarded as formed at the initial moment of time. 
Then one can write for G: 



ri tot (w)dwt 3 (-) 3 

-oo T 

+ r dw /V0 3 exp(-^+^^^(-) 3 lK,^^ 

Jw+ JO Vtotiw*) T $ 

The second term is negligible in comparison with the first one. Then 

Vtot (w)dwt 3 (-) 3 

-oo T 

and the behavior of supersaturation is determined. Then 



r] tot {w)dwt 3 {-) 3 /n c 

-oo T 
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The last step to do is to determine w*. For w* we have 

2 exp( 7o vtotM exp(i $ w 

or 

- =exp(- / — r exp(-- / r] to t{w)dwt (-) )dt) 

Having calculated the integral one can come to 

- = exp( — — / ?7 tot (w - ] 1 0.9) 

2 Vtot{w*) rioo® J-oo r 

The last equation together with an explicit expression for / given by the 
classical theory of nucleation is an ordinary algebraic equation. 

The number N^iw) of the droplets formed on the centers with activity 
w can be calculated as 

at f \ t \n t ( .I(w*,®,r} tot (w*)) 
Noo{w) = r)tot{w){l - exp(-exp(-u> + w*) — - 

The total number of droplets can be calculated by the integration over 
all w. 

Consider the opposite situation. The spectrum of activities participating 
in the nucleation process is so narrow the one can consider all centers as the 
similar ones. Then the system of equations will be the following 

i^,vtot) , r#) 

V(t) = »fe*exp( — I exp(-— -)dt) 

g{t) =m Vtot ) ['(t-t'f^eM- 1 ^)—^ 

JO T A 71^9 1] tot 

One can use iterations defined as r) i+ i = Q(gi), ft+i = P(gi,r)i) with 
initial approximations g = 0, i] = i] tot . Here Q, P define the r.h.s. of two 
previous equations. The second iteration gives already suitable results. Then 
one can get with a relatively high precision for the total number of droplets 

iV(oo) = [1 - exp[-0.9 (-p^a - ) 11 

The last relation solves the problem. 
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